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Recent data obtained in e + p collisions at HERA by the HI collaboration are presented: searches for new 
phenomena beyond the Standard Model, measurement of the inclusive neutral and charged current cross sections 
and of the proton structure function, determination of the gluon density in the proton and photon from jet 
productions, measurement and interpretation of event shapes in the current region of the Breit frame, test of the 
parton evolution at low x and of the diffractive production mechanism. 
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1. Introduction 

The electron-proton collider HERA offers the 
unique possibility to explore the complex struc- 
ture of the proton in a relatively clean envi- 
ronment. The large centre of mass energy of 
y/s w 300 GeV allows a perturbative QCD anal- 
ysis in a previously inaccessible domain. Com- 
pared to fixed target experiments HERA extends 
the kinematic reach of deep-inelastic scattering 
(DIS) by about of two orders of magnitude in the 
Bjorken variable x and the photon virtuality Q 2 . 
The understanding of the proton structure is in 
many aspects a challenge to theory and therefore 
interesting in itself. Furthermore, the proton is 
the particle which gives access to the highest pos- 
sible centre of mass energy in particle collisions. 
The knowledge of the (non-perturbative) parton 
density function in the proton and the theoreti- 
cal understanding how they evolve in the different 
kinematic regimes is therefore an important pre- 
requisite for any possible future discovery of new 
phenomena, e.g. in proton proton collisions at 
LHC. 

Within the Standard Model many aspects 
of the strong interactions remain unexplained. 
HERA has here the peerless opportunity to link 
QCD to the general analytic properties of the 
scattering amplitudes as e.g. exploited in the 
Regge phenomenology. Such theories have been 
very sucessful to describe elastic and diffractive 
scattering at high energies and the energy depen- 
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dences of the total hadronic cross section. First 
steps to understand the dynamics of the transfor- 
mation of a complex coloured system in a colour 
singlet probed by a virtual photon have been very 
sucessful. 

In the period 1994- 1997 HERA colliding E e = 
27.5 GeV positrons on E p — 820 GeV protons 
has delivered an integrated luminosity of about 
J C dt = 37 pb" 1 of useful data to HI. Most of 
the results in this report will be based on this data 
set. Since 1998 HERA has switched back to elec- 
trons and about J C dt = 15 pb _1 have been col- 
lected. First results based on J C dt = 5.6 pb _1 
from electron-proton collisions are already avail- 
able. This represents an increase in e~p data by 
a factor of 10. In addition the proton energy was 
pushed to E p = 920 GeV leading to a higher 
centre of mass energy and to a higher parton lu- 
minosity at fixed x. 

2. Search for New Phenomena at Highest 

Energy 

The high momentum transfers accessible at 
HERA allow the proton to be probed at very 
small distances down to 10 -16 cm via exchange 
of highly virtual gauge bosons. Particles with 
masses higher than the centre of mass energy 
could be exchanged in the ^-channel. Their inter- 
ference with standard DIS processes could be ex- 
perimentally observed as enhancement or deficit 
in the NC cross section at high Q 2 or at high 
mass. Moreover, new particles with masses up 
to 300 GeV could be produced as s-channel res- 
onances. HERA is here particularly sensitive to 



particles coupling to electrons and quarks in the 
proton. 

2.1. Inclusive Searches 
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Figure 1. M = \fsx distribution for the inclusive 
e + p — > e + X sample and for a optimised cut in y. 
Shown are HI data and the DIS expectation. 

In an analysis of the data taken in 1994 to 1996 



an excess of neutral current data e + p 
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at very high Q has been found. 12 events with 
Q 2 > 15000 GeV 2 have been observed in the 
data while only 4.71 ± 0.76 were expected by 
standard DIS. 7 of these events had an invariant 
mass M = ^/sx in a window of 25 GeV around 
M = 200 GeV and y > 0.4. In 1997 the total 
luminosity was doubled. Also in this data a ten- 
dency for more events at high Q 2 has been found. 
The magnitude of the excess is however not con- 
firmed. Above Q 2 > 15000 GeV 2 22 events 
are found and 14.7 ± 2.1 are expected. 8 events 
fall in the mass window at M = 200 GeV and 
3.01 ± 0.54 are expected. 

The mass distribution of all e + p data is shown 
in Fig. [|. The standard DIS background is known 
within ±5 — 10%. If the new particle has spin 
0, the angular distribution of the decay prod- 
uct is flat. The y = 1/2(1 + cos6>*) distribu- 



tion]^ will therefore be flat, as opposed to DIS 
where the cross section falls like 1/y 2 . The sig- 
nal to background ratio can therefore be opti- 
mised by a cut in y. The mass distribution for 
a ycut(M) optimised for each mass bin, is also 
shown in Fig. PH. From this mass spectrum upper 
limits on scalar and vector leptoquarks classified 
according to their spin, weak isospin and fermion 
number [ p| have been derived. If the strength 
of this new e + q coupling is the same as for the 
electromagnetic coupling, than masses up to 275 
GeV are excluded. If one allows for a branching 
LQ — ► eq of e.g. 10%, masses up to 255 GeV are 
ruled out by the data. For these small branch- 
ings this limit extents far beyond the domain 
presently covered in pp collisions. This has trig- 
gered intensive searches for other signatures. No 
signal has been found in the reactions e ± p — v \i p 
and e p — > r ± p. However, the derived limits on 
leptoquarks decaying in the second or third gener- 
ation are better as or competitive with low energy 
experiments in most cases [ |2J . 

2.2. Muon Events 

In an inclusive search for events with a high 
transverse momentum of the hadronic final state 
above 25 GeV, six events with high energetic 
isolated leptons (B£ > 10 GeV) have been ob- 
served. One event has an electron and five have a 
muon in the final state. In all events an imbalance 
of missing longitudinal and transverse momentum 
is measured suggesting the presence of an unde- 
tectable particle like a neutrino. The electron 
event and two muon events are found in the phase 
space where the most important background - the 
direct production of W bosons ep — ► WX - are 
expected. In the muon (electron) channel 0.8±0.2 
(2.4±0.5) events are expected and 5 (1) have been 
observed. 

The striking characteristics of the muon events 
is the large amount of transverse energy of the 
hadronic system (Py ) which is displayed together 
with lepton- neutrino transverse mass (Mj, ) in 
Fig. g. Overlayed is the expected distribution of 
direct production of W bosons. It would be how- 
ever interesting to know how QCD effects which 



2 0* is the lepton-parton angle in the centre of mass system 
of the decaying particle. 
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Figure 2. Transverse momentum of the hadronic 
final state as a function of the transverse mass 
of the lepton-neutrino system for events with iso- 
lated leptons. The crosses gives the la uncer- 
tainty on the measured parameters. The expected 
distribution for direct W production and muon 
production in 77 processes is overlayed. 



are presently not included in the background es- 
timation using the event generator EPVEC [ [|, 
could alter the P^f and Mlt v distribution. In 
the 1998/1999 e~p data no events with isolated 
leptons have been found. For a luminosity of 
JCdt = 5.1 pb- 1 0.14±0.04 (0.37±0.07) events 
are expected in the muon (electron) channel. 

3. Inclusive Neutral and Charged Current 
Single Differential Cross Sections 

3.1. Neutral Current Events 

The Standard Model neutral current cross sec- 
tion can be represented in the form: 



dV 



2-KOL Z 



dxdQ 2 xQ 4 



Y+F 2 (x,Q 2 )-y 2 F L (x,Q 2 ) 
T Y-xF 3 ( Xl Q 2 )] 



where a is the electromagnetic coupling constant 
and Y± = 1 ± (1 — y) 2 contains the helicity de- 
pendence of the electroweak interaction. 

The generalised proton structure function 
-^2(2, Q 2 ) can be decomposed in: 
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where F^™ describes the pure photon exchange, 
F2 the jZ interference and F% the pure ex- 
change of the Z boson. F 2 is sensitive to the sin- 
glet sum of the quark distributions (xq+xq). xF$ 
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Figure 3. Single differential neutral current DIS 
cross section as a function of Q 2 . 
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Figure 4. Single differential charged current DIS 
cross section as a function of Q 2 . 



is the parity violating proton structure function 
which is sensitive to the non-singlet difference of 
the quark densities (xq — xq) and is given by: 



F,= 



Q 2 



M 2 / 3 



1 z 



Q 2 



Q 2 



Ml 
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-F3 contributes with different sign to the cross 
section for the e~ and e + case. Below 
Q 2 < 1500 GeV 2 contributions from Z exchange 
are smaller than 1%, but at high Q 2 the inclusive 
cross section is sensitive to electroweak effects (at 
Q 2 = 5000 GeV 2 and x = 0.08: 6 Z -S 3 & 10%). 

The longitudinal structure function Fl con- 
tributes less than 5% at very high y and is negli- 
gible below y < 0.4. 

The single differential cross section is shown 
as a function of Q 2 in Fig. [| for both e + p and 
e~p scattering. The e~p data have been taken in 
1998/1999 and correspond to the / Cdt « 5 pb" 1 
analysed up to March 1999. 

The cross section falls like 1/Q 4 over 6 orders 
of magnitude and spans 2 orders of magnitude in 
Q 2 . For Q 2 > 3000 GeV 2 the e~p data are al- 
ways above the e + p data. This provides for the 
first time evidence for ^Z interference effects in 
e p scattering. The increase of the centre of mass 
energy from y/a = 300 to 320 GeV does not in- 
fluence this conclusion, since it leads only to a 
marginal increase of the cross section (see dashed 
line in Fig. 0). While in the e~p case good agree- 
ment is found with the Standard Model expec- 
tation, in the e + p case a slight overshoot is ob- 
served at the highest Q 2 . The theoretical uncer- 
tainty mainly introduced by the parton density 
functions is about 7% at the highest accessible 
Q 2 values. 

3.2. Charged Current Events 

Charged current events are characterised by 
their missing transverse momentum due to the 
undetected neutrino in the final state. Their cross 
section is given by: 
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where Gf is the Fermi constant and M prop is the 
mass of the space like charged exchanged boson 



which is in the SM the W boson. In leading order 
the helicity weighted parton density function have 
the simple form: 

*£?(*, Q 2 ) = (u + c) + (l- yf ■(d + s + b) 

$ e c£(x, Q 2 ) = (u + c) + (1 - yf -(d + s + b) 

The positrons (electrons) only interact with 
quarks or anti-quarks with negative (positive) 
electric charge. In the e~p (e + p) case mostly 
(anti-)quarks are involved in the interaction. 

Fig. [| shows the charged current cross section 
as a function of Q 2 . The e~p data are an order 
of magnitude above e + p data. This is due to the 
fact that in the e~p case the cross section is pro- 
portional to (u + c) while for e + p the coupling to 
valence quark is suppressed by (1 — y) (d + s). 
Moreover is the d quark density smaller than the 
u quark density in the proton. 

The charged current cross section falls over 
four orders of magnitudes in Q 2 . From the 
Q 2 dependence the propagator mass has been de- 
termined to be: 



M 



,„ „ P = 81.2 ± 3.3 (stat.) ± 4.3 (syst) GeV 



This result is compatible with the world average 
of the time-like W mass M w = 80.41 ± 0.1 GeV 

[1- 

At high Q 2 the charged current cross section 
is suppressed relative to the neutral current cross 
section due to the d/u ratio being less than unity. 
From this observation it has been determined that 
at X = 0.3 the u valence quark density is approx- 
imately 1.5 — 3 times larger than d valence quark 
density. A determination of u/d ratio at large 
x free from uncertainties due to nuclear binding 
effects will be possible in the future. 

These measurements represent a benchmark for 
the Standard Model of electroweak interactions. 
More e~p and e + p data at the highest possible 
Q 2 are needed to get the best sensitivity. HERA 
will increase the luminosity in the year 2000. 
High statistics data will then allow a better un- 
derstanding of the detector and decisive tests for 
deviation from the standard model are possible. 
This physics is just starting at HERA and a sur- 
prise is possible! 



Measurement of the Proton Structure 
Function 



Figure 5. The proton structure function F 2 as 
a function of Q 2 for bins in the Bjorken scaling 
variable x. Superimposed is a NLO QCD fit. 

A classic key measurement for the understand- 
ing of the proton structure is the precise deter- 
mination of F 2 by counting inclusively the lepton 
scattered off the proton. The kinematic domain 
spans about 5 orders of magnitude in Q 2 from 
1 < Q 2 < 30000 GeV 2 and covers the momen- 
tum range from the sea quark region at low x 
(x>2 • 10~ 5 ) to the valence quark region at large 
x (x < 0.65). 

An inclusive measurement has the advantage 
that it can be directly compared to QCD calcu- 



lations. However, small interesting effects on top 
of the dominating parton evolution are difficult 
to reveal. Therefore F 2 has to be precisely pined 
down to get a handle on the QCD evolution and 
to constrain the non-perturbative parton density 
functions. A summary of the presently available 
data is shown in Fig. [5| For x < 0.01 the sta- 
tistical error is about < 1% and the systematic 
accuracy has reached < 3—4% in most of the cov- 
ered range. The largest experimental systematic 
errors are introduced by the hadronic and electro- 
magnetic energy scale uncertainties. Compared to 
previous HI results [||, |], Q] the systematic uncer- 
tainties improved by nearly a factor of two. This 
precision can compete with the fixed target data. 
The measurement extents to energies of the 
scattered positron down to 3 GeV. This was 
only possible by using the backward electromag- 
netic fibre-lead calorimeter SPACAL (operated 
since 1995) and by using the backward silicon 
strip detector (operated since 1997) to measure 
tracks within 172° < Q d < 177°. It gives ac- 
cess to the kinematic region of low Q 2 and high 
y < 0.75 where the DIS cross section is sensitive 
to Fl. A new method has been introduced to 
improve the sensitivity to Fl by measuring the 
derivative: 

0F 2 _ ^ 2- 1 + dF L l_ 
Slog 2; Y? cHogx Y + 



da 



d log y 



at fixed Q 2 . Assuming that da/dlogy is a linear 
function of log y up to large y, the contribution of 
F 2 at high y has been extrapolated by a straight 
line fit for y < 0.2. The extracted Fl extents 
down to Q 2 —3 GeV at x rs 10~ 4 and is in 
good agreement with previous HI results [|| and 
consistent with Fl as calculated by NLO QCD. 

The NLO fit shown in Fig. g is performed 
in the MS renormalisation scheme using the 
NLO DGLAP evolution equations for three light 
flavours with the charm contribution added in 
the massive scheme according to a NLO calcu- 
lation of the boson-gluon fusion process. The in- 
put distributions of the valence (u v ,d v ) and sea 
(S = u = d — 2s) quarks and the gluon are pa- 
rameterised at a starting scale Q 2 = 2 GeV 2 us- 
ing 5 parameters which are determined in a fit to 
the data. The muon-proton data of BCDMS [ 



g] and the muon-deuteron data of NMC [ [l(| 
are included in the fit to constrain the high x 
behaviour of the parton distributions. Target 
mass corrections were applied for the fixed tar- 
get data. The strong coupling constant was set 
to a s (Mz) — 0.118 and the charm quark mass 
to m c =1.5 GeV. This fit gives a good descrip- 
tion of the data for all Q 2 and x. At low x the 
behaviour of F 2 is dominated by the gluon and 
the sea quark distributions. No need for a parton 
evolution different from NLO DGLAP is found. 

5. Jet production in DIS 

The high centre of mass energy of HERA leads 
to a large phase space for hadron production and 
to the possibility to observe clear jet structures 
in DIS. Jets are infrared and collinear safe ob- 
servables defined by an algorithm relating the un- 
observables quarks and gluons to the sprays of 
hadrons observed in the detector. Events with 
two jets in the central part of the detector can be 
produced in quark ((77 — > qg) or gluon (177 — > qq) 
induced hard processes. The leading order Feyn- 
man diagrams for dijet production are shown in 
Fig. J6[ Dijet cross sections are directly sensi- 
tive to the strong coupling constant a s and to 
the gluon densities in the proton. They can be 
calculated in pQCD to NLO in a s using numeri- 
cal methods implemented in several Monte Carlo 
programs. In a framePL where the virtual photon 
7 and the proton p collide head-on, it is obvious 
that two scales can characterise the hardness of 
the process: Q 2 as used in the QCD analysis of 
the inclusive DIS cross section or the transverse 
momenta of the two jets Et as used in hadron- 
hadron collisions. 

HI has reported the first fundamental under- 
standing of jet production in DIS at HERA en- 
ergies at the DIS 98 conference by demonstrat- 
ing that dijet cross sections for 10 < Q 2 < 
5000 GeV 2 and 8.5 < E T < 35 GeV can be 
described by NLO QCD. Jets were defined by 
the inclusive Kt algorithm [ O, [l^] in the Breit 
frame requiring two jets with Et > 5 GeV and 
Et,i + Et,2 > 17 GeV. Only jets lying well 
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Figure 6. Leading order feynman diagram for di- 
jet production in DIS at HERA. 

within the detector acceptance — 1 < ryjetjab < 2.5 
are considered. 

At this conference further progress has been 
reported. The double differential single inclusive 
jet cross section do-/dQ 2 dEx for 7 < Et < 50 
GeV has been made available. It is observed 
that for higher Q 2 the Et spectrum gets harder. 
Moreover dijet rates l/Ndn/dy 2 , where y 2 is the 
maximal value of the resolution parameter which 
allows to resolve two jets in an event, have been 



measured for 150 < Q < 



3 e.g. the hadronic centre of mass frame defined by 7+p : 
or the Breit frame defined by 7 + 2 x p = 0. 



5000 GeV 2 using 
the Durham Kt algorithm [ O with a scale of 
100 GeV 2 . Also the y 2 spectrum gets harder 
as Q 2 increases, i.e. harder jet structures are 
resolved. Both observables are well described 
by NLO QCD. Furthermore it has been demon- 
strated that the dependences of the internal jet 
structure on the Et and 77 is well reproduced by 
QCD models [@. 

The dijet cross sections d 2 a/dQ 2 d£ and 
d 2 a/dQ 2 dx have been used to extract the gluon 
density in the proton at a factorisation scale of 
fif = 200 GeV 2 which corresponds roughly to 
(E T ). £ is defined as£ = x (1+s/Q 2 ), where ^/§is 
the invariant mass of the dijet system. For Q 2 > 
200 GeV 55% of the dijet cross section is caused 
by gluon induced processes. In the fit a s (Mz) is 
assumed to be a s (M z ) = 0.119±0.005 [|§. This 
value has been mostly determined in processes 
which are independent of the proton structure. 
Inclusive DIS cross section at 200 < Q 2 < 650 
GeV 2 are simultaneously fitted. These data 
strongly constrain the quark densities, but de- 
pend only weakly on the gluon density. The re- 
sulting gluon density together with its error band 
is shown in Fig. |t] in the range 0.01 < x < 0.1. 
The error band is dominated by the uncertainty 
on a s , by the renormalisation scale dependence 
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Figure 7. The gluon density in the proton as a 
function of x determined in a NLO pQCD fit to 
dijet cross sections. 

and by the experimental absolute hadronic en- 
ergy scale. Good agreement with indirect deter- 
minations via scaling violations of F^ is found, 
but larger x- values are reached. Gluon determi- 
nations from global fits [ O, U7L n8[ tend to give 
slightly lower values, but are consistent within the 
errors. 

The precise knowledge of the gluon density at 
large x and of a s is an important ingredient to 
predict cross sections at the LHC [ Q . The re- 
gion 10~ 3 < x < 10 _1 is directly relevant for 
a Higgs with mass 100 — 500 GeV produced at 
\rj\ < 2. For instance, the cross section for the 
process gg — + Higgs gives within 20% different re- 
sults when different PDF parameterisations are 
used. Q 2 has to be extrapolated by DGLAP by 
3 orders of magnitudes from the HERA region 
to reach values relevant at LHC. Therefore the 
precise knowledge of the parton evolution is also 
a key ingredient in the cross section prediction. 
Also the cross section gg — > W which can be used 
as a luminosity monitor in proton-proton colli- 
sions varies within 5% when the most recent PDF 
of MRS and CTEQ are used. However, it is clear 



that comparing different PDF parameterisation is 
rather inadequate to estimate their uncertainty, 
since similar data are used in the fit and similar 
assumptions are made. Therefore, it is also im- 
portant to develop techniques to extract QCD pa- 
rameters together with an error estimate. Within 
one experiment HI has demonstrated that this is 
perfectly possible and more results are expected 
in the future. 

6. The Structure of the Real Photon 

7 „7 „7 




Figure 8. Feynman diagram for resolved processes 
in dijet production at HERA. 

In the quark parton model DIS is viewed as 
a highly virtual photon interacting with partons 
freely moving in the proton. This is a good ap- 
proximation when small distances are probed, i.e. 
in the limit where Q 2 is large. In this regime, the 
photon behaves like a point-like object, i.e. it di- 
rectly couples to quarks to produce the hard scat- 
tering. At low virtualities the photon dominantly 
fluctuates into vector mesons. However, the pho- 
ton may also fluctuate into a gq state with higher 
transverse energy without forming a bound state. 
In this case the photon acts as a source of strongly 
interacting partons. Such a 'resolved' process can 
be calculated within pQCD using the concept of 
a photon structure function. 

They can be experimentally investigated by de- 
manding jets or charged particles at high trans- 
verse energy in the final state. Such data are 
sensitiv to the quark and to the gluon content 
in the photon. Of particular interest is the be- 
haviour of the photon at low fractional energies 
of the parton participating in the hard scattering 
(x 7 ). 

The access to low x 7 is experimenally difficult. 
Since x 1 is given by x 7 = Yli2^T e~ n /2 _E 7 , 
low jet Et destroying a good correlation with the 
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Figure 9. -Dzjei cross section as a function o/x 7 . 

hard partons or large n values limited by the de- 
tector acceptance and the understanding of soft 
underlying processes have to be used. In a recent 
analysis of photoproduced events (Q 2 ~ GeV 2 ) 
jets are selected using a cone algorithm [ g(J with 
R = 0.7. After correction for energy steming 
from soft processes migrating into the jet cone 
("pedestal") Et > 6 GeV is required. Both jets 
have to lie well within the detector acceptance of 
—0.5 < rjj e t < 2.5 and their rapidity difference 
must not exceed |Ary| < 1. The dijet cross sec- 
tion is shown as a function of x 1 in Fig. 0. It is 
compared to the LO prediction of the PHO JET [ 
H || QCD model using the GRV-LO PDF [ || 
for the photon and the proton. While at i 7 ~ 1 
the direct contribution dominates, for x 1 < 1 
the resolved contribution is most important. At 
the lowest x 7 values gluon initiated processes are 
mainly responsible for the measured cross section. 
A leading order interpretation of the data can 
be made by using the concept of an effective PDF 
for which the dijet cross section can be written as [ 
|2l: 



V^ feS ( X 7' ^T> Q ) /, 

i,3 



j/p 
eff 



\ME, 






P /e (y,Q 2 ) 
v 

where the sum runs over all scattering matrix el- 
ements depending on the scattering cos 8* of the 
processes shown in Fig. H folded by the density of 
parton j in the proton and parton i in the pho- 



ton. The effective probability to find parton i in 
the photon can be written as: 
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The unfolded effective gluon distribution is shown 
in Fig. no as a function of x-y. The gluon density 
rises as x^ decreases. This rise is less pronounced 
than predicted by the LAC1 structure function [ 
M, but is in agreement with the GRV expecta- 
tion [ |23|]. The results improve a previous mea- 
and agree with earlier results us- 
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ing the complementary approach of single charged 
particles [ 071] . 
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Figure 10. Extracted leading order gluon density 
in the photon as a function of x 1 . 

7. Event Shapes in Current Region of 
Breit Frame 

Measurements of event shape variables F pro- 
vide information about perturbative and non- 
perturbative aspects of QCD. They allow to fit 
analytical expressions to data without referring to 
a fragmentation model by exploiting their char- 
acteristic power behaviour. Event shapes have 
been extensively studied in e + e~experiments at 
different center of mass energies [ E8| . 

Results from e + e~can be compared to DIS 
in the current hemisphere of the Breit frame. 
Thanks to the large kinematic range covered at 
HERA, the dependence of mean event shape val- 
ues on a hard scale can be studied in one exper- 
iment. Event shapes like thrust, jet broadening 
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Figure 11. Distribution of the C -parameter for 
various Q — y/Cp bins. 

and the jet mass have been investigated in the 
past [ M. In this workshop for the first time 
also the C-parameter and differential jet rates 
yi for the JADE and for the Kt jet algorithm 
have been presented. The C-parameter is defined 
as: C := 3(AiA2 + A2A3 + A3A1) where the A^ are 
the eigen values of the linearized momentum ten- 
sor of the final state particles. As an example, 
the distribution of the C-parameter in bins of the 
momentum transfer Q = \/Q 2 ranging from 7 to 
100 GeV is shown in Fig. |n|. While at high Q 
the data are described by a 0(a 2 s ) calculation [ 
P0[ , at low Q the shape does not agree with per- 
turbative QCD. The same conclusions hold for 
the mean value of the C-parameter shown as a 
function of Q in Fig. p~2. 

The mean event shapes can be expressed by a 
perturbative and a non-perturbative contribution 
of the form [|l|,|§: 

(F) = (F) pcrt + (F) pow 

(F) peTt = Cl , F (x)a s (Q) + c 2>F (x) a s 2 (Q) 

{F y™ = a ™!p:i n pQ [a o(Mi )_ as (Q)_ 



O 2 

(bln-^ + k 
Mi 



26 N 



l(Q)} 



where b = (11 C A - 2/)/12tt and k = 
[(67 — 37r 2 ) Ca — 10 /] /36-7T and ao is a free, but 
'universal', effective coupling parameter below an 
'infra-red' matching scale Aqcd <C /ii ^C Q. 

Such an ansatz is able to describe the data (see 
Fig. H3). The power corrections are large at low 
Q, but become less important with increasing Q. 
The parameters a s and ao can be simultaneously 
fitted to the data. It turns out that the analyt- 
ical form of the power correction is adequate to 
describe the data for all studied event shapes. A 
simple form (F) pow ~ 1/Q is not able to describe 
the data. 
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Figure 12. Mean value of the C-parameter as a 
function of Q = \j Q 2 ■ 

The question whether consistent fit results are 
obtained when using different event shape vari- 
ables is addressed in Fig. |l3|. The ellipse displays 
the 1<7 and 2cr contours of the statistical and the 
systematical error. The NLO dependence on the 
renormalisation scale is not included. The ex- 
perimental errors are highly correlated. The two 
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Figure 13. a s (Mz) versus arj(Mi) from the fit of 
four different event shapes variables. The ellipse 
gives the la and la contour of the statistical and 
the systematical error. 

parameters ao(p,) and a s (Mz) come out about 
in the same range like found in e + e~collisions. 
However, in particular the large spread of a s 
is not fully satisfactory. More precise data and 
more theoretical work is needed to continue this 
promising way to get an (analytical) understand- 
ing of hadronisation for specific variables. Open 
questions are the correct inclusion of higher or- 
ders to get a more consistent picture and the 
treatment of the uncertainty introduced by the 
parton density function to the fit result. 

8. QCD parton dynamics at low x 

At low x the simple picture of DIS as a process 
where a virtual photon interacts instantaneously 
with a point-like parton freely moving in the pro- 
ton has to be modified. The phase space for gluon 
emission (W 2 « Q 2 /x) between the photon and 
the proton becomes so large that many partons 
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Figure 14. Diagram of a DIS events at low x. 

can be radiated before interacting with the pho- 
ton. Such an interaction is illustrated in Fig. |IJ. 
In particular in the central rapidity region the 
detection of hard particles can discriminate be- 
tween different evolution schemes of the parton 
cascade [|33j]. 

The transitions g — > gg and g — > qq at each 
point in the ladder can be approximated by the 
DGLAP equations [ [}4J. By resumming the 
(a s lnQ 2 ) terms they predict the Q 2 evolution 
of a parton known to be point-like at some given 
scale Qq to the region where the interaction with 
the photon takes place. To derive the DGLAP 
equations a strongly ordered configuration in the 
parton virtualities along the ladder has to be as- 
sumed. This leads to a suppression of the avail- 
able phase space for gluon radiation towards the 
proton. 

When the (a s lnl/:r) terms become large, 
they have to be taken into account, e.g. by the 
resummation accomplished by the BFKL equa- 
tions [|3j|. In a physical gauge, these terms corre- 
spond to an n-rung ladder diagram in which gluon 
emissions are ordered in longitudinal momentum. 
The strong ordering of the transverse momenta 
is replaced by a diffusion pattern as one proceeds 
along the gluon chain. The BFKL equations de- 
scribe how a parton in the proton is dressed by a 
cloud of gluons at low x localised in a fixed trans- 
verse spatial region of the proton. 

The density of hard partons in the central 
rapidity region can be experimentally explored 
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Figure 15. ir° cross section as a function ofx. 

by identified particles or jets in the forward re- 
gion of the HI detector (see Fig. Qj). The 
finely granulated HI LAr calorimeter with its 
3.5x3.5 cm calorimeter cells and a 4 fold segmen- 
tation offers the opportunity to identify high ener- 
getic 7T° mesons by exploiting the typical proper- 
ties of electromagnetic showers. The high particle 
density in the forward region for 5° < 6i a b < 25° 
requires a good understanding of the detector and 
is the main experimental challenge. The inclu- 
sive cross section as a function of x for 7r° mesons 
with a F T * > 2.5 GeV and x n = E n /E p > 0.01 
in the centre of mass system for three Q 2 regions 
ranging from 2.0 < Q 2 < 70 GeV 2 is shown 
in Fig. [15| The cross section rises by a factor 
8 towards low x for each Q 2 . In the lowest Q 2 
region values down to x = 5 • 10 -5 are reached. 
A weak dependence of the cross section on Q 2 is 
found and for each Q 2 the x spectrum is similar. 
A Monte Carlo study indicates that a ir° with 
B£ rs 2 GeV originates on average from a parton 
with P T * « 5 GeV. 

The LEPTO Monte Carlo [ |f§ model based on 
the leading order matrix element (0(a s )) eq — ► 



eqg or eg — ► eqq incorporating higher order emis- 
sions from the proton side in a leading logarithm 
approximation close to the DGLAP equations is 
not able to describe the data. In the highest 
Q 2 bin LEPTO is below the data by a factor of 
2. This difference increases towards lower Q 2 . 
The RAPGAP model [ || includes in addition 
LO resolved photon processes. It is in much bet- 
ter agreement with the data. However, at the 
lowest x and Q 2 also this model fails to describe 
the data. An analytical calculation based on the 
leading order BFKL equation [ p7| agrees with the 
data. In this calculation the parton calculation is 
related to the measurement using a fragmentation 
function. 

Whether this measurement can be taken as a 
proof that BFKL effects are needed to describe 
HERA data remains an open question. It is not 
fully transparent what approximations have been 
included in the calculations in view of the large 
NLO correction to the BFKL equation. Moreover 
it is not a priori clear if the fragmentation func- 
tions are valid in the high particle density en- 
vironment at HERA. Furthermore, it has been 
shown that very similar measurements of jets in 
the forward region [ |3£| £39[ can be described by 
a 0{a 2 ) calculation of direct and resolved pro- 
cesses [ |4t|. Since the photon splitting term is 
the dominant contribution of the resolved part, 
this calculation suggest that a 0(oP s ) direct calcu- 
lation would be enough to describe the hadronic 
final state data at low x. In this case the parton 
ladder is not long enough such that resummation 
effects to all order of the In (1/x) terms are needed 
in the HERA regime. 

The strong rise of the inclusive 7r° cross section 
reminds the rise of F^ towards low x which has 
been before the advent of HERA predicted by the 
BFKL equations. We know today that F 2 can be 
well described by the NLO DGLAPequations, but 
is also compatible with BFKL. It is interesting to 
note that the ratio of the forward ir° cross sec- 
tion to the inclusive cross section is constant over 
the full x range. At 2 < Q 2 < 4.5 GeV 2 ap- 
proximately 0.2% of the DIS events contain a for- 
ward 7T° with B£ > 2.5 GeV. At 15 < Q 2 < 70 
GeV 2 this ratio rises to 0.5%. While for an inclu- 
sive measurement the DGLAP evolution is ade- 
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Figure 16. Diagram for a colourless exchange in 
DIS. 

quate, for observables where the parton evolution 
is probed at a specific r\ region the implemen- 
tation of the DGLAP equation in Monte Carlo 
models are too restrictive. Only models allowing 
for a large phase space of parton emission are in 
agreement with the data. 

The constant (with respect to the inclusive 
cross section) probability to emit hard particles in 
the forward region can be juxtaposed to the prob- 
ability to emit no particle at all. It is known since 
1993 that the ratio of events exhibiting a large ra- 
pidity gap is approximately 10% for 5 < Q 2 < 50 
GeV 2 and2-10~ 4 < x < 2-10~ 2 . Does this mean 
that there is a connection between the production 
mechanism of rapidity gap events and the parton 
dynamics at low x? 

9. Diffraction 

A new era in understanding hadronic interac- 
tions was opened by the observation that a sur- 
prisingly large fraction of DIS events (w 10%) 
exhibited a rapidity region free of hadronic activ- 
ity between the particles emerged from the hard 
scattering and the proton. Since in normal DIS 
large rapidity gaps are exponentially suppressed, 
these events were soon analysed in terms of an ex- 
change of a colourless object probed by the virtual 
photon. If no quantum numbers are exchanged 
such interactions are called "diffractive" and the 
colourless object is called "pomeron" (P) 

The diagram for the production of diffractive 



events is depicted in Fig. [16]. The hadronic final 
state can be split up in two distinct systems X 
and Y which are separated by the largest rapid- 
ity gap in the event. Usually events are selected 
by requiring an absence of activity in the forward 
part of the detector such that Y has a small in- 
variant mass (My < 1.6 GeV) and the squared 
momentum transfer t between the incoming pro- 
ton and Y is small (|i| < 1 GeV 2 ). The following 
kinematic variables can be defined to describe the 
diffractive production mechanism: 



&■■ 



Q 2 + M 2 



Q 2 + Ml 
Q 2 + W 2 



Q 2 



M 2 



Q 2 



M x 



x p measures the fraction of the proton momen- 
tum transfered to the P and is the momen- 
tum fraction of the P momentum carried by the 
quark coupling to the virtual photon. If the sys- 
tem Mx is further resolved in a subsystem M12, 
the variable z^ gives the momentum fraction of 
the struck parton. If M\ 2 < M x , i.e. z w < 1, 
not the entire energy from the colourless exchange 
takes part in the hard scattering. 

Usually the diffractive interaction is described 
in terms of Regge phenomenology and the kine- 
matic dependences are parameterised as: 



da 



dtalM x 



oc (W 



2\ 2a ip(*)- 2 e M 



where a w (t) is the effective leading Regge trajec- 
tory. For soft diffractive processes, a w (t) takes 
the universal form a^it) ~ 1.08 + 0.2 t. When a 
hard scale is present, 0^(0) (t-slope) is expected 
to increase (decrease). 

The "soft" pomeron ansatz successfully de- 
scribes the total and elastic hadron-hadron and 
hadron-photon cross section as well as the en- 
ergy dependence of the cross section of light vec- 
tor mesons (Mx = My). However, a single soft 
pomeron does not describe all diffractive data 
measured at HERA. In presence of a hard scale 
like Q 2 or a larger vector meson mass the cross 
section steeply increases with W lp . A new type 
of dynamical pomeron may begin to play a role 
whose structure can be tested in DIS. HERA 
allows diffractive interactions to be analysed in 
terms of parton dynamics and offers the possibil- 
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ity to study the transition between perturbatively 
calculable and incalculable strong interactions. 

9.1. Inclusive Diffractive measurements 

The inclusive diffractive data are consistent 
with a partonic interpretation of the pomeron for 
which the cross section factorises into the parton 
density of the pomeron and a pomeron flux fac- 
tor describing the probability to find a pomeron 
in the proton [ Bl| . The parton densities can be 
determined by a QCD analysis [ E3] of the inclu- 
sive diffractive cross section: 



d a ev ^ e xY 



' ep- 



47ra 



dQ 2 d/3dx n 



PQ 



Td-y+\)F? i3 \ X]P ,Q 2 ,P). 



At fixed Xjp the Q 2 dependence of the diffrac- 
tive structure function F 2 ^ '(x F ,Q 2 , j3) is rather 
flat. The j3 dependence exhibits a logarithmic 
rise characteristic for scaling violations. When 
fixing the quark and gluon distribution at a start- 
ing scale of Qq — 3 GeV 2 and evolving them 
to larger Q 2 using the DGLAP equations, it is 
found that most of the partonic content of the 
pomeron is carried by hard gluons. At Q 2 = 4.5 

GeV 2 the pomeron consist to 90% of gluon and 
even at Q 2 = 75 GeV 2 the gluon fraction is still 
80%. Using these fits performed for Q 2 < 100 

GeV 2 also the new HI data reaching Q 2 up to 
800 GeV 2 can be described. However, in an 
analysis of very low Q 2 data ranging down to 1 

GeV 2 the logarithmic (3 dependence flattens off 
around Q 2 «3-4 GeV 2 . 

9.2. Diffractive Dijet Production 

The partonic structure of diffractive interac- 
tions can be further tested in the hadronic final 
state. High transverse momentum jet production 
is directly sensitive to the gluon content. For 
7.5 < Q 2 < 80 GeV 2 and 0.1 < y < 0.7 ex- 
actly two jets are selected with Exjet > 5 GeV. 
Jets are defined by the cone algorithm. Exjet is 
measured relative to 7* axis in rest frame of the 
system X. The dijet cross section as a function 
of z w is shown in Fig. [IT]. Models in which the 
pomeron is dominated by gluons (fit 1 and 2 in 
[ |42|) describe the data well. Models where the 
pomeron consist only of quarks at the starting 
scale undershoot the data by a factor of 5. 
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Figure 17. Diffractive jet cross section as a func- 
tion of Zjp . 

An alternative model sees the diffractive inter- 
action as qq fluctuation of the photon (decom- 
posed in Fock states) coupling via two gluons to 
the proton [ |||, Q in a colour singlet configura- 
tion (2-gluon model). At the largest z w values, 
where the full momentum of the X system is car- 
ried by the two jets and the model is expected to 
be valid, it is somewhat below the data, but still 
compatible with them given the large systematic 
error. 

In conclusion, dijet data support the gluon 
dominated partonic picture of the pomeron. 

10. Diffractive Charm Production 

The role of gluons in the diffractive produc- 
tion mechanism can be directly assessed by tag- 
ging charmed quarks in the hadronic final state, 
since they predominately originate from Boson- 
gluon processes. Charm quarks are identified by 
reconstructing D*^ mesons in the classical de- 
cay chain D* -> D°tt+ 1ow -> {K-n+)^ ow . The 
branching ratio of this process is only 2.6%. In 
the kinematical region of 2 < Q 2 < 100 GeV 2 
and 0.05 < y < 0.7 45 events have been se- 
lected in a data sample of about 21 pb" 1 taken 
in 1995 — 1997. D* mesons were required to have 
Pt > 2 GeV and \n\ < 1.5. The total cross sec- 
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Figure 18. Diffractive charm cross section as a 

function of z w . The pomeron model has been di- 
vided by a factor of 3, the soft colour model by a 
factor of 2. 
tion is measured to be: 

a ms (ep^ (D*X)Y) = 154±40(stoi)±35(syst)pb. 

In Fig. |l^ the visible cross section a V i S (ep — > 
(D*X)Y) is shown as a lunction ol z w . As for 
the diffractive jet data the cross section rises to- 
wards low z w . This means that in most of the 
events a part of the hadronic final state produced 
in the diffractive exchange is not associated with 
the hard subprocess. 

The shape of the z w distribution is well repro- 
duced by the gluon dominated partonic pomeron 
picture. This model remarkably reproduces the 
rise of the cross section towards low z w . How- 
ever, the normalisation is off by a factor of 3 ! 
Hence, much less charm quarks are produced in 
the hadronic final state than predicted by the 
parton density function extracted from the inclu- 
sive diffractive measurements and applied to the 
pomeron picture. 

In the region where the hadronic system X 
dominantly consists of the cc pair alone (i.e. 
z w ~ 1) the 2-gluon model gives a fair descrip- 
tion of the data. As in the case of the diffractive 
dijets it is not able to describe the the low z w 
region. The soft colour interaction model where 
the colour structure of a normal DIS events is 
modified by a soft interaction leading to a colour 
singlet in the final state [ E5[ can not describe 



the overall normalisation nor the shape of the z w 
distribution. 

In conclusion, the diffractive charm signal is 
more pronounced at low z w values. The data 
therefore indicate that not the entire colourless 
exchange couples to cc system. The behaviour 
is well reproduced by the "pomeron" model, but 
the overall measured cross section is three times 
lower than expected. This is the first time that 
this picture which has been so far supported by 
all inclusive and hadronic final state data is put 
into question. More data are needed to clarify the 
situation. 

11. Vector Mesons 

Elastic: 




j/t or p 
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Figure 19. Diagram for vector meson production. 

HERA allows to study the simplest photon- 
hadron diffractive reaction ^*p — ► VM p in a 
new range of "f*p centre of mass energies up to 
Wy P < 300 GeV. It is possible to investigate the 
transition from the soft regime of photoproduc- 
tion of light vector mesons at low momentum 
transfers to the pQCD regime in the presence of 
three hard scales: Q 2 , Mx = My and t. The tra- 
ditional way of modelling vector meson produc- 
tion is the vector meson dominance model where 
the photon fluctuates into a vector meson which 
then elastically scatters off the proton. In such 
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a picture only a weak energy dependence of the 
cross section is expected [ Q (see Fig. 19 upper 
left). This works well for light vector mesons at 
low Q 2 and t. For heavier vector meson like the 
J/ip the charm mass provides a hard scale and 
pQCD calculations are possible. An example for 
such a pQCD model where the J/ip couples via 
two gluons to the proton is shown in Fig. [L9] (up- 
per right). 



failure of this picture together with the sucess of 
QCD inspired models gives additional confidence 
in the potential of QCD to describe processes 
where the scales involved are surprisingly small. 
Further studies will gain a better understanding 
of the transition from the region governed by soft 
interactions to the regime where pQCD turns on. 
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Figure 20. The elastic J/ip cross section as a 
function of W lp . 

The measured elastic J/ip cross section as a 
function of W lp is displayed in Fig. |2(] for < 
Q 2 < 33 GeV 2 . A wide W 1P range from 
50 < Wy P < 285 GeV is covered. A much steeper 
rise (Wyp' 77±0 ' 18 ) than predicted by soft models 
(Wyp' 32 ) is observed in the data. However, the 
data are well described by a QCD inspired calcu- 
lation [ |4?| . The arrow indicates the uncertainty 
introduced by the charm mass (m c = 1.4 — 1.5 
GeV). Since the gluon density enters squared in 
the calculation this measurement is sensitive to 
the gluon distribution in the proton. The cross 
section decreases as 1/(Q 2 + M 2 ,,) , where n = 
2.38±0.11, as Q 2 increases. The data are not pre- 
cise enough to see a further steepening of the rise 
with increasing Q 2 . In conclusion, VDM holds 
unless there is a hard scale in the problem. The 
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Figure 21. Differential cross section for inelas- 
tic J/ip production as a function of Q 2 , P^, Z, 
Wyp and the rapidity. 

Inelastic vector meson production is charac- 
terised by an inelasticity < z = (P<p ■p)/(q-p) < 
0.9 and a relatively large mass Mx- Two models 
for the production meachnism are shown in the 
lower part of Fig. [19|. In a short distance process 
a cc state is produced in a boson-gluon fusion. 
The cc system is therefore a colour octet. The 
long range transition into the J/ip has therefore 
to involve the emission of an additional gluon. 
The failure of this approach to describe the pro- 
duction of quarkonia at large Pt, in particular 
at TEVATRON, led to the inclusion of additional 
contributions where the cc system transforms into 
an colour singlet in a long range interaction with- 
out emitting a gluon [ Bq] . The long range ma- 
trix elements have been fitted to the pp data. It 
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is therefore interesting to see, if these fits can be 
sucessfully applied to HERA data. 

The differential cross section for inelastic 
J/tp production defined by Mx > 10 GeV is 
shown in Fig. KlL The colour octet model is over- 
layed as dashed the singlet model as solid line. 
Both predictions [ Q fail to describe the shape of 
the data. The colour singlet is a factor 2 — 3 below 
the data, the sum of the singlet and octet model is 
too large. The reason for this failure is unknown. 
It seems that an overall adjustment of transition 
matrix clement is necessary. The failure to de- 
scribe the shape also calls for an relative adjust- 
ment of individudal contributions. Recently it 
has been reported that the inclusion of higher or- 
der radiation (in a Monte Carlo approach) lowers 
the colour octet matrix elements needed to de- 
scribe the TEVATRON data [|(J. Similar effects 
might also be important at HERA. 

12. Conclusions 

In the past years HERA has been operated 
very successfully and provided useful e p data 
sets. From the recent 1998/1999 e~p running pe- 
riod neutral and charged current DIS cross sec- 
tion have been measured. The observed differ- 
ence in the e~p and e + p data provides clear evi- 
dence for electroweak interference effects. HERA 
starts to probe the electroweak sector of the Stan- 
dard Model. So far, no striking deviations have 
been seen, although for Q 2 > 15000 GeV 2 more 
events than expected are found in the e + p data. 
Future high precision data will reveal, if this is 
only a statistical fluctuation or if it is a first hint 
for a new interaction on top of QCD. The HERA 
collider will be upgraded in the year 2000 to de- 
liver an integrated luminosity of 150 pb _1 per 
year. This will allow to accumulate J Cdt sa 
1 fb _1 until the year 2006. 

Precise data on the proton structure func- 
tion constrain the parton density functions in 
the proton. This inclusive measurement is con- 
sistent with the parton evolution according to 
the DGLAP equations. The gluon density has 
been directly measured from observables based 
on the hadronic final state like dijet and charm 
production. One of the biggest uncertainties in 



these measurement is the unknown relation be- 
tween the partonic final state and the measur- 
able hadrons. First steps to develop analytical 
calculations for specific observables based on the 
hadronic final state in the current hemisphere of 
the Breit system are very encouraging. Further 
theoretical and experimental work is needed here. 

New tests of forward particle production put 
the DGLAP evolution scheme into question. 
More partonic activity seems to be required. 
Moreover it has been observed that the proba- 
bility to emit hard partons in the forward region 
is constant with respect to the inclusive cross sec- 
tion. This is very similar to the behaviour of event 
where no particles are emitted in a rapidity re- 
gion. 

So far HERA data have supported a consistent 
picture of the production mechanism of rapid- 
ity gap events. New results on diffractive charm 
production seem to be in disagreement with the 
pomeron picture. More data are needed to clarify 
this situation. 

Precise data on vector meson production are 
now available. In the elastic channel QCD mod- 
els are able to describe the data. However, in 
the inelastic channel no satisfying calculation is 
available. 
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